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1 Introduction 

In this paper, we consider harmonic maps from weighted graphs into metric 
spaces. This can be considered as a generalization of geodesic lines in Rieman- 
nian manifolds. Our considerations are based on the simple observation that 
a geodesic considered as a map j : [0, 1] — ► N from the unit interval to the 
Riemannian manifold N and parametrized proportionally to arclength is char- 
acterized by the property that for all sufficiently close < a < b < 1, j( s ^-) is 
the unique midpoint of 7(a) and 7(6), that is, 

= argimn qEN (d 2 (j(a),q) + d 2 (j(b),q)). (1) 

This leads us to represent a geodesic as a string of points in N each of which 
is the midpoint of its two neighbors. At the same time, this allows us flexible 
refinements, that is, we can insert additional points in the string as midpoints 
of consecutive ones already present. For that, it is useful to also consider the 
following slight generalization of {]]) 

7(ta + (1 - t)b) = argmin 9eAr (^ 2 ( 7 (a), q) + (1 - t)d 2 ( 7 (6), g)) (2) 

(0 < t < 1). 

A midpoint is a center of gravity of two points. In a Riemannian manifold, such 
centers of gravity exist locally uniquely, that is, when the points whose center is 
to be constructed are sufficiently close. Globally, uniqueness need not be true. 
Therefore, we may need to localize in the image. 

It is then clear how to conceptualize a harmonic map from a weighted graph 
into N. We simply require that the images of the nodes of the graph are 
appropriately weighted centers of gravity of their neighbors. Here, in order to 
localize in the image, we might need to refine the graph by subdividing edges. 
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Harmonic maps from graphs into compact Riemannian manifolds were studucd 
in [8]. Our approach, however, naturally leads to a generalization to metric 
space targets that locally admit such unique centers of gravity. This class of 
metric spaces includes the important class of Alexandrov spaces with upper cur- 
vature bounds, see [T] as a systematic reference. 

Thus, in this paper, we show the existence of harmonic maps from weighted 
graphs into such metric spaces, shortly called harmonic nets. (When we have 
a space with non-positive curvature in the sense of Alexandrov or Busemann, 
this result is contained in a general existence result for harmonic maps of the 
first author, see [HE]-) The proof is not difficult. It is based on the iterative 
replacement of image points by the centers of gravity of the images of their 
neighbors, following the strategy described in [7], together with suitable adap- 
tive refinements to keep the constructions local. Here it is important that the 
domain, that is, our graph, can be treated as a one-dimensional object. While 
two dimensions represent a border line case, in higher dimensions, general con- 
structions of harmonic maps are only possible when the target space possesses 
non-positive curvature. The reason that is that the energy functional we are 
employing is quadratic, and therefore the scaling behavior is different in dimen- 
sions 1, 2, and greater than 2. Therefore, the essential features of our scheme 
are local uniqueness and the scaling property of our functiona. 
Our constructions possess certain similarities with some schemes employed in 
numerical analysis, like the standard difference scheme for the numerical solu- 
tion of the Laplace equation or adaptive refinements in multigrid methods. A 
key conceptual feature of our approach is that we systematically exploit local 
uniqueness of solutions and that we need to make explicit only the images of a 
discrete set of points because then all other images are implicitly determined by 
that local uniqueness. Therefore, as in good numerical schemes, we never have 
to work out or store more information than we actually need. 

2 Geometric concepts 

We let (N, d) be a complete metric space. For abbreviation, we usually simply 
write N in place of (N,d), the metric d(., .) being implicitly understood. We 
say that N admits refinements if for any p, q E N, there exists some m 6 N 
with 

d(m,p) = d(m, q) = -d(p, q). (3) 
We also call such an m a midpoint of p and q. 

Definition 2.1. Suppose that N admits refinements. We define the radius 
r(N) of unique refinement as the largest r £ [0, oo] with the property that for 
any two p,q € N with d(p, q) < 2r, their refinement (midpoint) m = m(p, q) is 
unique. 
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More generally, we say that q £ N is a center of gravity of the finitely many 
points pi, ...,p n £ N with weights w\, ...,w n (> 0) if 

n 

q = argmin p ^ w j ^(P'Pj)- ( 4 ) 

3=1 

We define the convexity radius c(N) as the largest c £ [0, oo] with the property 
that whenever p,pi, ...,p n are points in N with d(p,pi) < c for i = 1, ...,n, then 
the center of gravity of pi, ...,p„ (with any positive weights) is unique. 
Since a midpoint of two points is their center of gravity when they are given 
equal weights, we have obviously < c(N) < r(N). 

We let r be a finite weighted graph with vertex set I and edge set E where 
each e £ E has a weight w(e) > 0. We say that two vertices i, j are neighbors if 
they are connected by an edge. Thus, for our purposes, a graph is a discrete set 
I together with a symmetric neighborhood relation ~ and (symmetric) weights 
w(i,j) — w(e) for any two neighboring vertices i, j (i.e. i ~ j) connected by the 
edge e. 



We define the refinement r r of T as the graph with vertex set I U E with i £ I 
and e £ E being neighbors when i £ e in T; there are no other pairs of neighbors 
in r r . The weight of such a pair of neighbors is w(i, e) = \/2w(e) where the 
latter is the weight of the edge e in T. The edge set of T r then is obvious. We 
can then also define successive refinements of T r . 

A map / from r to N assigns to every i £ I some point p = f(i) in N. 
We define the energy of such an / : T — ► N as 

E(f) = j2m), (5) 

iei 

with E i (f)= ™ 2 (i,j)dV(i)J(j))- (6) 

jel; i~j 

In particular, for i ~ j, 

d\m,m) < ^-fi(f)- (7) 

We say that the map / is harmonic if for all i £ I, f(i) is a center of gravity 
of the points /(j), j ~ i, with weights Wj — w 2 (i,j). 
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3 Characterization by angles in tangent cones 



The above concepts of refinement and center of gravity find their natural place 
in the context of Alexandrov's metric spaces. For a systematic development of 
this theory that we shall use in this section, see pQ . 

These spaces enjoy particular properties when their (Alexandrov) curvature is 
bounded from above. It is part of the definition of such a space of curvature 
bounded from above that any two sufficiently close points can be joined by a 
shortest geodesic which then is in fact unique and depends continuously on these 
endpoints. (We may also parametrize it by arclength - and call it an arclength 
geodesic - if convenient.) Some of the general notions in the theory, however, 
do not need the assumption of an upper curvature bound. That assumption 
then is rather employed to derive geometric properties of the objects defined in 
the theory. 

An important concept here is the tangent cone of a metric space at a point. 
Let (N, d) be a metric space, 71 , 72 : [0 , e] — ► (N, d) be arclength geodesies 
emanating from a point P 6 N . 

Consider points Q € 71, R € 72 different from P. An (upper) angle 0(71,72) 
between 71 and 72 is defined as 

— d 2 (P 7 Q)+d 2 (P 7 R)~d 2 (Q,R) 
cos0( 7l ,7 2 ) :=hm Q ^P 2d(P ,fl) d (p, Q) ■ 

We have the following characterization of the angle between 71 and 72 (see [2] , 
[II. 1-II.3] ) : let (N, d) be a metric space whose curvature is bounded from above 
by a constant K > 0. Then 

co S %i,7 2 ) = to mjm^ki^km , 

s^O s 

provided in case K > that e is less than the diameter of the comparison model 
space of constant curvature K. 

A geodesic curve 7 starting at a point P 6 N has a direction if 0(7, 7) = and 
two curves have the same direction if the angle between them is equal to zero. 
This is an equivalence relation on the space of curves starting from the same 
point P £ N and the completion of the set of equivalence classes (endowed with 
the distance induced by the angle) is called the space of directions f2 p (N) of N 
at the point P. 

The tangent cone TpN of (N, d) at a point P 6 N is the cone over the space of 

directions, namely T P N = (Sl P (N) x R+) / (Q P (N) x {0}). 

We will denote a tangent element by [7,3c], where 7 € £lp(N), x > and 

elements [7, 0] are identified with the origin O v of TpN . 

The distance d in N induces on TpN a distance function dp defined by 



dp{[li,xi] , [72,^2]) = 



x\ +x\ - 2xix 2 cos6(^ ll 72) if #(71,72) <tt 

Xl + X 2 if #(7l,72) > 7T 
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For those [7, x] for which we can find a unique geodesic from P with direction 7 
that can be extended up to distance x, we define that point as the exponential 
image of [7, x]. The inverse of this exponential map, the projection map from 
the subset of N where it is defined to the tangent cone TpN is denoted by np. 
In the case of (simply connected, complete) nonpositively curved metric spaces, 
it is well known that the map 7rp is defined everywhere, distance non-increasing 
and distance preserving in the radial direction (see |10j). 
The following important result has been proved by Nikolaev [S]: 

Lemma 3.1. Let (N,d) a metric space of curvature < K, with K > 0. Then 
the tangent cone at a point of N is a space of non-positive curvature in the sense 
of Alexandrov. 

Let P, Q, R be points in N and Q s = (1 — s)P + sQ the point on a distance 
realizing geodesic joining P and Q with distance s.d(P, Q) from P. 
We have the following Taylor expansions: 

d 2 (Q s ,R) = d 2 (P, R) - 2sd(P, R) cos P (Q. R) + a(s) , with lim ^ = 

s^O s 

d P (ji P {Q s ),i:p{R)) = d 2 {P,R)-2sd{P,R) cos 9 P (Q,R)+b(s) , with lim — = , 

s^O S 

where 9p(Q, R) denotes the angle subtended by Q and R at P. 

We recall the following result concerning harmonic maps (see [3l 110]): 

Proposition 3.1. Let f : T — ► (N,d) be an harmonic map, then: 

(i) : Vie/, 7r/(i)(/(i)) minimizes the function w 2 {h j)^(t)(- > n fli)f(J)) 

in Tfu\N. 

(ii) : Vie/, E^ 2 (U) (V, 7r /(i) /(j)) < , V V e T /(i) iV, 

where ( , ) denotes the inner product defined on Tfu\N by: 
([71, xi] , [72,22]) = xi x 2 cos 0(71,72) 

(m) : Vie/, i/ie barycenter in //(,) 0/ £/ie points (7T/(i)/(j))j'~i wif/l 
weights { w ,j?$ coincides with the origin Oi := TTf(i)f(i) ofTf^N, where 

w(i) = x> 2 («, i)- 

The inequality (ii) in the above proposition will be interpreted as the critical 
condition for harmonic nets. 
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4 Refining maps 



If N admits refinements, we can construct a refinement f r : T r — > N of a map 
/ : r — > N by assigning to every edge e connecting i and j in T a midpoint of 
/(i) and /(j). We observe that for each i G T, we have 

Wr) = ^(/) (8) 
where on the left hand side, i is considered as an element of r r . Also, 

£>(/r) = X> e (/r) = il5(/r) (9) 

by symmetry, where we consider the is and es as vertices of T r . In particular, 
we have from ©, © 

Lemma 4.1. 

£(/r) = £(/)■ (10) 

-V / : r — ► -/V is harmonic then so is its refinement f r . 

The converse holds when distances between images are sufficiently small, that 
is, when midpoints between the images of neighbors are unique. 

From ([7|) and (|10[1 , we conclude that by performing sufficiently many successive 
refinements, we may assume that all distances between the images of any two 
neighboring vertices are smaller than some prescribed e > 0, for example smaller 
than r(N) or c(N) when that quantity is positive. 

5 Homotopy classes 

For the present purposes, we write r(f) and r(T) instead of f r and T r because 
we wish to consider the refinement as an operation that can be iterated. For 
example, r 2 (r) = (r r ) r is obtained as the refinement of r r . A refinable map 
/ : r — * N then is considered as a collection of iteratively refined maps r n (f) : 
r n (T) -► N for n e N. 

Assume now that the refinement radius r(N) > 0. We say that two maps 
/l) f% '■ T — > N are geodesically close if for every i 6 T, d{f\{i), < 2r(N), 

that is, the images of i under fx and f% have a unique midpoint. A refinement 
of the pair /i,/2 then is defined to be the triple /i,/i,2>/2 where /i^ is the 
midpoint map of f± and fa, that is, for every i S T, /i,2(i) is the midpoint of 
/i(i) and / 2 (i). 

We say that two refinable maps f,g:T—>N are geometrically homotopic if 
there exists a sequence /o = /, /i, /2, /a = .9 for some i e N of refinable 
maps such that for any n £ N and any 1 < j < A, the maps r n (fj-i) and 
r n (fj) are geodesically close. This sequence can again be refined by putting in 
midpoint maps between consecutive sequence elements. 

Geometric homotopy is an equivalence relation, and the equivalence classes are 
called geometric homotopy classes of refinable maps from T to N. 
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6 Construction of harmonic nets 



We assume that N admits centers of gravity. By subdividing suitable edges of 
r as above, we may assume that T is bipartite, that is, its vertex set is a disjoint 
union / = I\ U I2 such that all the neighbors of any point in one of those subsets 
are contained in the other one. On the space C = C(T, TV) of maps / : T — ► N, 
we define maps p a : C — ► C, a = 1, 2 with p a {f) being the map obtained from 
/ : r — > N by replacing the image of every f(i) for i 6 I a by a center of gravity 
of the f(j) for j ~ i. As long as the centers of gravity are not unique, we need 
to make choices here, but in the situation where c(N) > 0, we can assume that 
r has been sufficiently refined (depending on an upper bound E for the energy 
of /) so that the images f{j) of the neighbors of any i E T possess a unique 
center of gravity. (This follows from ([7]) and the fact that the edge weights get 
multiplied by a factor of V2, that is, become larger, under each refinement.) In 
that case, the maps p a (f) are unambiguously defined for all / with E(f) < E. 
p a decreases (or, more precisely, does not increase) the energy density Ei(f) 
for all points in I a , but not necessarily for those in the complement of I a . 
Nevertheless, since by symmetry £ ieJl = E ie / 2 ^(/) = \ E if) ( scc 

([9])), we have 

Lemma 6.1. 

E(p a (f)) < E(f) (11) 

for all f . 

Moreover 
Lemma 6.2. 

E{pMf))) = E{f) (12) 

if and only if f is harmonic. 

Theorem 6.1. Let (N,d) be a compact metric space that admits centers of 
gravity. Let T be a finite weighted graph. Then, for any map f : T — > N , the 
iterations f n := (P2P1)™/ contain a subsequence converging to a harmonic map. 

Proof. Since N is compact, we can find some sequence v(n) of positive integers 
going to infinity for which f u i n ){i) converges to some point fo(i) G N for every 
vertex i of the finite graph T. We have 

fu(*+i) = {P2piY [n) /„(„) for some p{n) G N. (13) 

Since the metric d behaves continuously under convergence (since it defines the 
topology of N) , we have 

E(f ) = lim £(/„(„)). (14) 

n — >oo 

But then also 

lim £(/,,(„+!)) = ]imE((p 2Pl )^f v{n) ) 

< lim E{f v r n \), by Lemma I5TT1 because of /x(n) > 1 
= E(f ). 
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Thus, equality has to hold throughout. Moreover, pip\f v {n) converges to 
and so 

E{p 2 pih) = limE((p 2Pl r^ +1 U in) ) as before 

= E(fo) from the preceding observation. 

Lemma 16.21 then implies that fg is harmonic. □ 

The assumption of the theorem that the space N admits centers of gravity is 
satisfied when N has an upper courvature bound. For k £ K, we denote by 
Dk the diameter of the n-dimensional, complete, simply connected model space 
with constant sectional curvature k. We then have the following result from 
Alexandrov's theory (see [I]): 

Lemma 6.3. Let X be an Alexandrov space with curvature bounded above 
by k. Then, for every x G X , there exists a positive number R x £ (0, fy] such 
that the closed metric ball centered at x and with radius R x , B(x,R x ), is a 
convex subset in X . 

Remark: When N has non-positive curvature in the sense of Alexandrov or 
Busemann, our theorem is contained in a general theorem of the first author 
(see [1]) and when TV is a compact Riemannian manifold, it follows from the 
fact that any homotopy class contains at least one harmonic map (see [5]). 
Since distances of neighboring image points are controlled by the energy of a 
map, see ([7]), we see that if the refinement radius r(N) is positive, we may control 
the geometric homotopy class by assuming an energy bound and sufficiently 
refining the graph T. 
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